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Abstract 

We calculate the two-loop QCD mixing of the dimension-five flavor-changing operators that 
arise in the Standard Model after integrating out the W, Z and Higgs bosons and the top- 
quark, i.e. the mixing among gluonic and photonic magnetic moment operators. Our calculation 
completes the two-loop anomalous dimension matrix of operators that govern low energy flavor- 
changing processes. It is an important ingredient of the next-to- leading calculation of the 
B — > Xgj decay rate. 
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1 Introduction 



Flavor changing processes at energies much below the electroweak scale are most conveniently 
described with help of the effective hamiltonian [|TJ 

^/~EW,M (i) 

i 

where Oj are dimension-five and -six operators built out of light (u, d, s, c and b) quarks, leptons, 
photons and gluons, and Cj(/x) are their Wilson coefficients. The operators 0{ that remain after 
applying the equations of motion (EOM) are dimension-six four-fermion operators and the 
dimension-five magnetic moment operators: 

7 = eqa^F^q' (2) 
O g = gqa^G%T a q'. (3) 

Here F^ u and G a denote the photonic and gluonic field strength tensors, respectively, and T a 
are the color group generators. The quarks q and q' have different flavor and usually have 
specific chiralities but, as explained below, this is irrelevant in our calculation. The QED and 
QCD coupling constants are denoted by e and g, respectively. 

The factorization of short- and long-distance QCD effects is achieved |l| by evolving the 
coefficients Cj(//) from the electroweak scale down to the low energy scale of the process under 
consideration. The anomalous dimension matrix that governs this evolution is known up to 
two loops for the QCD mixing of the phenomenologically important 4-quark operators into 
themselves 0, quark- lepton operators £| and the magnetic moment operators 0,[|,|6|]. The 
only remaining two-loop QCD mixing among all these operators is the two-loop mixing of the 
magnetic moment operators (§) and (H) which we calculate in the present paper. 

Our results can have many phenomenological applications, as they affect all possible flavor- 
changing processes involving two quarks and a photon or a gluon. The two-loop mixing we 
calculate can also be directly applied to flavor- conserving magnetic moment operators (see 
end of section ^). However, the most important application of our results is their use in the 
calculation of the complete next-to-leading logarithmic QCD corrections to the B — > X s ~f 
decay rate. The critical phenomenological importance of calculating the next-to-leading QCD 
corrections to the B — ► X s ^ rate has been stressed in the previous paper 0. We will not 
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repeat this discussion here. Instead, we will put more emphasis on the details of the present 
calculation. 

The paper is organized as follows: In the next section, we recalculate the one-loop mixing 
of 7 and O g . We recover also the nonphysical one-loop counterterms that will be needed for 
cancellation of sub divergences in the further two-loop computation. Section 3 presents certain 
details of the two-loop calculation. Our final results are given in the end of this section. In 
the last section, we rewrite our results in the notation common for the analyses of the b — > 57 
decay and show the size of the next-to-leading contributions we have calculated. 

2 The one-loop calculation 

Our starting point is the lagrangian density 

£ = £qcdxqed + X ^C 7 7 + CgOg + ^2 CkOkj ■ (4) 

Here Cqcdxqed is the full QCD x QED lagrangian for f masslessQ flavors, including the 
gauge fixing and ghost terms. X denotes some normalization constant that does not get QCD- 
renormalized. The operators Ok are additional dimension-five operators required as countert- 
erms for off-shell Green's functions with insertions of 7 and O g . According to the standard 
textbook by Collins ||, these operators either vanish by EOM or are BRS-exact (i.e. are BRS 
variations of some other operators), and are therefore nonphysical. In the following, we will 
assume nothing about these operators but just recover them from the explicit one-loop calcu- 
lation. 

We have performed our calculation off-shell in the background field || version of the 
Feynman-'t Hooft gauge. We have calculated UV-divergent parts of one-particle-irreducible 
(1PI) one- loop diagrams with single 7 or O g insertions. All such diagrams with nonnegative 
degree of divergence correspond to the following transitions: q' — > q, q' — > qj, q' — > qQ, q' — ► qA, 
q' — > g77, q' q^Q, q' ~^ QlA, q' — > qQQ, q' — > qQA and q' — > qAA, where Q and A denote 
the quantum and the background gluon fields, respectively. The diagrams corresponding to the 
first three transitions are shown in fig. |1|. 

Since QCD and QED interactions preserve chirality, it is obvious that the mixing we calculate 

is independent of whether the operators 7 and O g are taken as they stand in eqs. (|j) and 
introducing fermion masses would not affect the mixing among 7 and O g calculated here. 
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Figure 1: One-loop diagrams for the q' — > q, q' — > g7 and g' — > gQ transitions. The squares 
denote insertions of 7 or O g . 

(0) or specific (opposite) chiralities for the quarks q and g' are chosen, as it happens in the 
phenomenologically important case of the b — > s'j decay. In the following, we will take these 
operators as they stand in eqs. (0) and (|3|). 

We find that divergences in all the possible one-loop 1PI Green's functions with single 7 
or O g insertions can be removed by counterterms proportional to 7 and O g themselves, and 
to the following two operators: 

Odd = 0H (5) 
Oqd = igq(-0 + Q0)q', (6) 

where is the usual QCD x QED covariant derivative. 

Assembling the four operators into the vector (0 7 , O g , Odd, Oqd) T , the corresponding one- 
loop anomalous dimension matrix reads 
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where N denotes the number of colors, Cf = (N 2 — 1)/2N and Q q is the electric charge of 
the g-quark (which must be identical to the one of q'). The mixing of 7 and O g among 



themselves, as well as the mixing of g into Odd are in agreement with what has been found 
previously []5|JT0f| . Note that in refs. [|5],|l(J and in most other papers, the quark mass is introduced 



in the normalization of 7 and O g . So 7h ere = 7there + 7ml > where 7 m = —-^-^GOf + 0{g i ). 

The mixing of O g into Oqd has not been considered so far because it arises from one- loop 
diagrams with external quantum gluon legs. We need to consider such diagrams as they are 
divergent subdiagrams in the further two-loop calculation. 

One may ask why only the very operators Odd and Oqd arise as nonphysical counterterms. 
We know || that all such counterterms either have to vanish by the EOM or be BRS-exact. 
Both Odd and Oqd vanish by the EOM for the quarks. In our case, there is not "enough 
dimension" for the EOM for the gauge bosons (such EOM has dimension 3 and the two quark 
fields have also dimension 3). So any EOM- vanishing operator must be proportional to the 
EOM for a quark. Two quark fields and a covariant derivative together have dimension 4. The 
remaining one unit of dimension can be provided either by another covariant derivative (then 
the only possibility is Odd) or by the quantum gluon field (as in Oqd)- The background gluon 
field cannot show up because our operators must be invariant under the gauge transformation 
given in eqs. (2.14) and (2.15) of ref . |gfl 
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8AI = -f abc u b Al - -d^ a (8) 

6Q a = _ f a b c ubQ c (g) 



where u a denotes the gauge transformation parameter. The two above transformations combine 
to the usual gauge transformation for the full gluon field + Q®. 

The final observation is that our initial operators 7 and O g are odd under charge con- 
jugation combined with exchange of q and q 1 . This is why the two terms in Oqd come in a 
combination odd under this symmetry. 

BRS-exact operators cannot arise as counterterms to 7 and O g . The arguments for this 
are the following: The BRS variation raises both dimension and ghost number by one unit. 
So any BRS-exact operator that could mix with 1 or O g should be a BRS variation of a 
dimension-four operator with ghost number —1. The only possibility for the latter operator is 
fj a qT a q' . The BRS variation of such an operator contains a term with two fermions and two 

ghosts. On the other hand, such a term cannot arise as a counterterm for the magnetic moment 
2 In our conventions, the sign of g is opposite to that in ref. M. 



operator because the corresponding diagram in fig. is convergent .[] 




Figure 2: Convergent diagram for the q' — > qr\r\ transition. 

As follows from the above considerations, one can show that 7 , O g , Odd and Oqd are 
the only possible counterterms for 7 and O g without calculating a single Feynman diagram. 
Then the matrix in eq. (0) can be found just from the diagrams depicted in fig. [I] with only 
quantum gluon on external lines. The calculation is then identical to the one in the usual 
Feynman-'t Hooft gauge (without background field). 

Since our calculation was fully computerized, it did not require much effort to check that 
indeed all the divergences in the diagrams with single 7 and O g insertions (also diagrams with 
two external gauge bosons) can be renormalized by the same counterterms proportional to 7 , 
O g , Odd and Oqd- This was one of the cross-checks in our calculation. 

The last remark in order is that the operator Oqd would arise also in the non-background- 
field calculation. Then Q in Oqd would stand for the full gluon field. Consequently, this 
operator would not be gauge-invariant (not even BRS-invariant). However, this is not a problem 
because EOM-vanishing operators are nonphysical anyway, i.e. their physical matrix elements 
vanish. They do not mix into physical operators, either ||. This is why the unknown entries 
in the lower right corner of the matrix in eq. (f7|) are completely irrelevant. 

One may then wonder why we have decided to recover of the mixing of O g into Odd and 
Oqd at all. The reason for this is the method we have applied for calculating the two- loop 
diagrams in the next section. 

3 The two-loop calculation 

At the two-loop level we need to calculate the q' —> q'j and q' — > qA Green's functions with single 
7 or O g insertions. Their divergent parts give the two-loop mixing among these operators. As 

a by-product, we obtain the two-loop mixing of O g into Odd- The latter mixing must come out 

3 At first glance, this diagram has zero degree of divergence (remember that the operator vertex has dimension 
1). But one power of momentum in the power counting is only the momentum of the outgoing external ghost 
— just because of the structure of the gluon-ghost vertex. So the effective degree of divergence is —1. This 
argument applies to all orders. 
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the same both from photonic and gluonic Green's functions. This is a nontrivial cross-check of 
the calculation. 

The main problem in the two-loop calculation is the large number of diagrams to be con- 
sidered. For instance, there are 167 two-loop diagrams for the self-mixing of O g . Therefore, 
a simple systematic algorithm is necessary to make the calculation fully computerized. With 
such an algorithm it is not really important whether one has to calculate 100 or 1000 two-loop 
diagrams — it is only a matter of computer time which in our particular case was around two 
days (on a SUN ELC workstation). 

The input to our programs were the Feynman rules following from the Lagrangian in eq. @ . 
The program FeynArts |TTJ was used for constructing the relevant topologies (see fig. ^), drawing 
the diagrams and creating the analytic expressions for the corresponding amplitudes. 





(11) 



(2) 



(7) 



(35) 



(18) 




(8) (3) (1) (10) (7) 

Figure 3: Topologies of the Feynman diagrams for the two-loop self-mixing of O g . The numbers 
in brackets are the numbers of diagrams corresponding to each topology. 



Next, a single algorithm (around ten pages of Mathematica [12] codeR) was applied to each 



of the diagrams. The potentially UV-divergent expressions were expanded to first order in 
external momenta and then the tensor integrals were reduced to scalar ones. So all the integrals 
we needed to calculate were scalar two-loop integrals independent of external momenta. 

There were no physical masses in our calculation. But before the expansion in the external 
momenta was made, an auxiliary mass was introduced in all the propagator denominators. Oth- 
erwise expanding in external momenta would create spurious infrared divergences that would be 
indistinguishable from the UV-divergences within the framework of dimensional regularization. 



4 We had two programs created independently by each of us. 
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These infrared divergences would not cancel in the sum of all the diagrams^] 

The auxiliary mass was the same for all the propagator denominators. So all the two-loop 
scalar integrals we needed to calculate were of the form 

dPq d D r 1 



(10) 



{2n) D {2n) D [q 2 - m 2 } a [r 2 - m 2 ] b [{q + r) 2 - m 2 ] c 

where a, b and c are integers and D is the dimensionality of spacetime. A change of variables 
(q — » q + r, r —>■ —r) shows that I a bc is symmetric under permutations of a, b and c. If any of 
these integers is nonpositive, the integral is easily reduced to a product of one-loop integrals. 
For a, b and c all positive one has 

,2' 



-^211 
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+ convergent terms 
+ convergent terms 



+ convergent terms, for 



(11) 
(12) 

n>3. (13) 



(4vr) 4 e (n - l)(n- 2)m 2 (™- 2 ) 

where e = 4 — When a, b and c are positive and two of them are greater than 1, then J a j, c 
is convergent. This fact and the above three equations are derived in the appendix. 

4])-(13) together with similar expressions for one- loop scalar integrals reduce the 



Eqs. 



calculation of the divergent part of any integral I a bc in eq. (|1(]) to simple algebra. The further 
steps of the calculation are of course purely algebraic, too: The results of the scalar integrals 
were combined back into tensor integrals and then contracted with the Dirac matrices, partly 
using Tracer |]T3. Performing the Dirac algebra was straightforward in our case, as we had 



only a single open fermion line and no 75 in our calculation. So no "evanescent operators" |14 
had to be introduced. 

Before presenting our final results, we should give arguments why introducing an auxiliary 
mass in all the propagator denominators and then expanding in external momenta leads us to 
the same results for the divergent part of the sum of all the diagrams^] as we would get without 
introducing this auxiliary mass and without expanding. 

The key point is the well-known fact that once all the sub divergences are removed in a 

mass independent renormalization scheme, the pole part part of a dimensionally regularized 

5 The simplest way to convince oneself about this is to consider an example with only one relevant diagram, 
e.g. one-loop wave-function renormalization in the scalar massless A^ 3 theory. 

6 i.e. the sum of the two-loop diagrams and the one-loop counterterm diagrams which remove subdivergences 



Feynman diagram is polynomial in masses and external momenta. This fact can be shown 
on the basis of given Feynman rules alone, and no field-theoretical arguments have to be used. 
The proof M is based on Weinberg's theorem |15|] and on the property, that differentiation with 
respect to external momenta decreases the degree of divergence of a Feynman diagram with no 
sub divergences. Adding auxiliary masses in propagator denominators does not affect this proof 
at all. So, after introducing these masses, the pole part of a Feynman diagram with subtracted 
sub divergences still will be a polynomial in external momenta and all the masses (including 
the auxiliary ones). Therefore, setting the auxiliary masses to zero in the end, we recover the 
result for the pole part of the diagram with no auxiliary masses. However, the latter is true 
only for the pole part and only after subtraction of sub divergences. 

The one-loop diagrams with counterterm vertices that cancel the subdivergences in two-loop 
diagrams have to be calculated with exactly the same Feynman rules as the two-loop diagrams. 
The same refers to the calculation of the one-loop renormalization constants themselves. So we 
also need to use massive propagators e.g. in all the diagrams for the one-loop renormalization 
of the two-point Green's function for the gluon. This necessarily leads to introduction of a 
"gluon mass" that gets renormalized at the one-loop level. It is known [TS] that a theory 
with a bare mass for a vector boson has no physical sense.f] But so long as the would-be 
gauge-fixing term is present in the Lagrangian, it is renormalizable by power counting [[17]], 
i.e. the Green's functions (even if physically meaningless) can be rendered finite with help of 



local counterterms. Moreover, in a mass-independent renormalization scheme (such as the MS 
scheme we use), these counterterms will be exactly the same as in the massless case (except for 
the renormalization of the vector boson mass itself). So the Feynman rules for the physically 
meaningless theory with a massive gluon can be used for recovering renormalization constants 
in true QCD. 

Similar arguments can be used to justify that we could introduce an auxiliary mass in all 
fermion propagator denominators without introducing it also in the numerators. Needless to 
say, this also affects the calculation of the "gluon mass" counterterm. 

Introducing auxiliary masses does not affect the results for the one-loop mixing summarized 
in eq. (0). Neither it affects the usual QCD renormalization constants (except for renormaliza- 
tion of mass terms). However, it affects the role the EOM- vanishing operators Odd and Oqd 

7 Unless this boson couples only to conserved currents. 
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play in our calculation. With the modified Feynman rules, the usual arguments which show 
that on-shell matrix elements of these operators vanish no longer hold. These operators remain 
nonphysical in the end, when we set the auxiliary masses to zero. But they are important at 
the intermediate stages of our calculation, i.e. the one-loop counterterm diagrams with coun- 
terterms proportional to Odd and Oqd do affect the mixing among the physical operators 7 
and O g . 

A remark concerning expansion in external momenta is now in order. As the two-loop di- 
agrams for q' — > g7 and q' — > qA transitions have degree of divergence equal to +1, the pole 
parts of our diagrams are linear in external momenta. This is why we expand in external mo- 
menta only to the first order before performing the loop integrations. The expansion in external 
momenta can be viewed as an exact splitting of propagators into parts that are polynomial in 
external momenta and parts that contribute to integrands with lower degree of divergence, as 
in the following example: 

\ - 1 | - 2 gP-P 2 1 (14) 

(q + p) 2 — m 2 q 2 — m 2 q 2 — m 2 (q + p) 2 — m 2 

where p denotes the external momentum and q is the loop momentum. Performing such an 
operation appropriately many times (in a diagram with subtracted sub divergences) one can 
split the integral into a convergent part and a part that is a polynomial in external momenta. 
This shows that for the pole part of a diagram with subtracted sub divergences, the expansion 
in external momenta is allowed before the loop integration is made. 

Our final results were subject to several cross-checks. First, the double pole part of the 
sum of all the two-loop diagrams had to be twice smaller and have the opposite sign than 
the double-pole part of the sum of all the one-loop diagrams with counterterm vertices. Only 
then the sum of all the diagrams together can be local (i.e. polynomial in external momenta). 
Second, the pole parts needed to have the structure of the gauge- invariant operators 7 , O g 
and Odd (the operator Oqd could not appear as we have not calculated two- loop diagrams 
with an external quantum gluon). This led to relations between coefficients at different Lorentz 
structures (usually two relations for four independent structures). Next, the mixing of O g into 
Odd had to come out the same from both q' — > q^ and q' — > qA diagrams. This is also a 
manifestation of gauge invariance. It is interesting to note that the latter requirement was 
fulfilled only after proper inclusion of the "gluon mass" one-loop counterterms. Finally, the 
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two-loop self mixing of 1 which is given by a relatively small number of diagrams has been 
calculated also "by hand" without introducing auxiliary masses and expanding in the external 
momenta. 

The final result of our two- loop calculation is the explicit expression for the 0(g A ) part of 
the matrix 7 in eq. (0). The anomalous dimension matrix now reads 
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(16) 
(17) 
(18) 

(19) 



Let us now briefly come to the case of two identical quark flavors, q = q' . Then we have to 
take into account several additional diagrams of the type depicted in fig. |j. All these diagrams 
contain a trace over an odd number of Dirac matrices (remember we are working with massless 
quarks) and so give no contribution. Thus, the results (p!5|)-([T9|) can be directly used for 
flavor-conserving magnetic moment operators as well.P| 



Figure 4: Examples of additional diagrams for flavor-conserving operators. 



3 We would like to thank Matthias Jamin for bringing this fact to our attention. 
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4 The b — > 57 case 

Having in mind the potential application of our results to the calculation of the next-to-leading 
logarithmic corrections to the B — ► X s 7 decay rate, we now rewrite the 2x2 anomalous 
dimension matrix of the magnetic moment operators once again in the conventional notation 
used in the analyses of this process. If we normalize the two magnetic moment operators as in 
ref. 

Qi = ^m b s a a^{l + l5 )b a F^ (20) 
Q 8 = £- 2 m b s a a^(l + l5 )(TT^G; u (21) 

then their 2x2 anomalous dimension matrix can be obtained from the corresponding entries 
in eqs. (p)|)-(|T6|) just by subtracting 7 m from the diagonal terms. At the two- loop level, the 



anomalous dimension of the mass reads [18[ with our sign convention 



7~ = -r^C, + J^C, (-'IN - 2D, + ) + 0(f). (22) 



Consequently, the anomalous dimension matrix for Q 7 and Q$ is 

.2 4 
167T 2 1 (167T 2 ^ 



7(i7) = 7 (0) t|^ + 7 (1) 7T#^ + 0(A (23) 



where 

7 (0) = i .... ! = \ i : i (2u 





and, in the MS scheme, 



7 (i) 



(7/ (f^iV - 16C, - f /) 
Q.Cf{&N-32C f -ff) -f -^ + ^ 2 + f|-f/iV 



4688 
27 

2192 4063 



27 



(25) 

81 27 / 

Unfortunately, knowing the above anomalous dimension matrix is not enough to write the 
full next-to-leading RGE for the coefficients of Qj and Q$. As explained in detail in ref. [^], 
the three-loop mixing of the four-quark operators into the magnetic moment operators gives 
contributions of the same order as the two-loop self-mixing of the latter operators. Therefore, 
our calculation is an important step towards the next-to-leading calculation of the B — > X s ~f 

11 



rate, but by no means completes this calculation. We hope that the methods described in the 
present paper can provide helpful tools in calculating the desired three-loop mixing. 

The numerical value of the next-to-leading corrections we have calculated can be easily 
determined from the matrix in eq. (Pq ) and the analysis presented in ref. 0. The equations 
(4)-(9) of that paper summarize the leading-order expressions for the b — > 37 decay rate which 
is proportional to the square of the quantity denoted there by C^ e ^ ' . Inserting the matrices 
from eqs. (p4]) and (|25|) above into the equations (21)-(28) of ref. j7|, we find the next-to-leading 



contribution that adds to C?^ in the MS -scheme: 



eff _ a s (M w ] 
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Z297664 _ x 7164416 _jl 256868 u 6698884 ie\ ™ /]f , 

+ 77 23 ri 23 -| 7723 7723 (M w ) 
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(26) 



with r/, Cf\M w ) and c£\M w ) defined as in eqs. (5), (8) and (9) of ref. 0. For a s {M z ) = 0.12, 
m t = 174 GeV and /x = 5 GeV, the above contribution causes a decrease of the predicted 
branching ratio by around 5%. However, no phenomenological conclusions can be drawn be- 
fore all the other next-to-leading corrections are calculated. Only then the total correction is 
renormalization-scheme independent. 

The leading order expression for the coefficient of the gluonic operator is given in the last 
equation of ref. [0. The next-to- leading contribution from the two- loop self- mixing of Q$ reads 



ACf" M = (,* - f*) Cf\M w ) (27) 



(m the MS -scheme) 

yeff(iA _ a s {M w ) 64217 
4tt 9522 

which in the Standard Model case gives numerically around 2.4% decrease of the coefficient, 
i.e. a similar change as in the photonic case. Even when the full next-to-leading calculation 
is completed, its effect on the coefficient of the gluonic operator is not expected to have any 
measurable physical consequences. So we should treat eq. (p7[) only as a by-product of the 
calculation performed for the photonic operator Qj. 

In the end, we would like to point out once again that our two-loop anomalous dimension 
matrix is insensitive to the scheme used for 75 in dimensional regularization. Two-loop mix- 
ing usually depends on the scheme used for 75, and the scheme dependence of the two-loop 
anomalous dimension matrix cancels in physical quantities with the scheme dependence of the 
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one-loop matrix elements In the case of Q7-Q8 mixing, we have completely eliminated 

7 5 from the calculation by using chirality conservation in QCD and QED interactions. This 
was equivalent to a trivial application of the recent Pisa scheme for 75 ||. However, exactly 
the same anomalous dimension matrix would be obtained in any of the other commonly used 
schemes HV, NDR and DRED (see Buras and Weisz in ref. || for a summary of these schemes). 
So the situation here is much simpler than in the case of four-quark operators. 

To conclude: We have calculated the two-loop QCD mixing of the magnetic moment opera- 
tors which are the only EOM-nonvanishing dimension-five flavor-changing operators that arise 
in the Standard Model after integrating out the W, Z and Higgs bosons and the top-quark. 
Our calculation is an important contribution to the phenomenologically desired next-to-leading 
QCD corrections to the b — > S7 decay. This was the last two-loop QCD mixing to be calcu- 
lated for the effective hamiltonian describing flavor changing processes at low energies in the 
Standard Model and many of its extensions. 

Acknowledgements 

We would like to thank Andrzej Buras, Matthias Jamin, Manfred Lindner, Uli Nierste, Hubert 
Simma and Peter Weisz for helpful discussions. This research has been supported by the 
German Bundesministerium fiir Forschung und Techno logie under contract 06 TM 732, by the 
CEC Science project SC1-CT91-0729 and by the Polish Committee for Scientific Research. 

Appendix 

In order to derive the equations (|11D-(|T3"D one starts from the integral 

J 211 (M, m lt m 2 ) = J ^ _ _ mw _ _L _ _ _ _ _ . (28) 

Differentiating the integrand with respect to mi or m 2 we find an integral that is convergent 
by Weinberg's theorem [|15j (i.e. all the subdiagrams of the corresponding scalar diagram have 
negative degree of divergence). Consequently, the divergent part of J 2 n(M, mi, m 2 ) is inde- 
pendent of m 1 and m 2 . We can therefore replace J 2 n(M, mi, m 2 ) by J 2 n(M, 0,0) and easily 
calculate the latter with help of the Feynman parametrization 

J 2 n(M, mi, m 2 ) = J 2 n(M, 0, 0) + convergent terms = 
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2 1/ , M 2 \] . . 

— H — 1 — — 2 In — — + convergent terms (29) 
(47r) 4 |_e 2 e \ 47r / J 

By integrating and/or differentiating the above result with respect to the masses appropriately 

many times and then setting all the masses equal, we obtain the equations (P1)-(|i~3|). This is 

also how one learns that I a t, c is convergent when two of its indices are greater than 1 and all of 

them are positive. 



The short derivation presented in this appendix is motivated by the appendix of ref. ITS 
where also the finite parts of the integrals are given. However, we disagree with the global sign 
of the explicit expression for (M 1 M\Mi\M 2 ) given there. 
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